B s -> KW°KW° DECAYS: 
THE GOLDEN CHANNELS FOR NEW PHYSICS SEARCHES 
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We point out that time-dependent CP asymmetries in B g — ► X*°K*° decays probe the presence 
of new physics in b — > s transitions with an unprecedented theoretical accuracy. We show that, 
contrary to the case of Bd — > <j>Ks, it is possible to obtain a model-independent prediction for 
the coefficient S(B S — > K*°K*°) in the Standard Model. We give an estimate of the experimental 
precision achievable with the next generation of B physics experiments. We also discuss how this 
approach can be extended to the case of B s — > K*°K° , B s — > K*°K° and B B — > K°K° decays and 
the different experimental challenges for these channels. 



The measurement of CP asymmetries in flavour chang- 
ing neutral current processes represents a crucial test of 
the Standard Model (SM). In particular, time-dependent 
CP asymmetries in b — » s penguin-dominated modes 
are considered among the most sensitive probes of New 
Physics (NP) 0. Measuring these asymmetries is one of 
the highlights of the B-factory physics program 0, [H, H| • 
In this context, the study of Bd — > <j>Ks has been consid- 
ered for a long time the golden mode for NP searches in 
nonleptonic B decays, since it is a pure penguin [5J. In- 
deed, writing the amplitude in terms of renormalization 
group invariant (RGI) parameters, defined in ref. [f|, one 
obtains: 



A(B d <f>K°) = v; b v ts P - V* b V us p 



GIM 



(1) 



where P contains penguin contractions of charmed 
current-current operators together with the matrix ele- 
ments of b — ► s penguin operators, while p GIM repre- 
sents the GIM-suppressed difference of penguin contrac- 
tions of current-current operators containing charm and 
up quarks respectively. 

Neglecting the contribution of P GIM on the basis of 
plausible dynamical arguments, the Bd — > 4>K$ decay is 
mediated by a single amplitude, so that under this as- 
sumption no direct CP violation can be produced and 
the time-dependent CP asymmetry probes the phase 
2/3 of Bd — Bd mixing. In terms of the coefficients S 
and C of sine and cosine terms in the time-dependent 
CP asymmetry, this means that in the SM one expects 
S(B d -> cf>K°) = sin 2/3 and C(B d -> <j>K°) = 0. The 
theoretical error associated to this prediction is related 
to the ratio of P GIM /P. 

The average of currently available experimental mea- 
surements by BaBar Q and Belle [3] gives S(Bd — > 
cj)K ) = 0.39 ± 0.18 and C(B d -> (/>K ) = 0.01 ± 0.13 7]. 
Even though the experimental errors are still large, it is 
interesting to observe that the value of S(Bd — ► 4>K°) de- 
viates from the world average sin 2/3 = 0.675 ± 0.026 Q- 



If confirmed in the future with a smaller error, this mea- 
surement might provide a hint of NP in B decays. On 
the other hand, this interpretation should consider the 
theoretical unccrtanty introduced by neglecting p GIM . 
Even though several model-dependent approaches, based 
on factorization or on flavour symmetries, have been pro- 
posed in literature [8j, |9[ , a model- independent evaluation 
of the error is not available yet. Considering that the next 
generation of B physics experiments (loj is expected to 
reduce the error on S(Bd — > <f>Ks) down to a few per- 
cent, the lack of a model independent evaluation of the 
theoretical error is a strong limitation for a complete and 
meaningful test of the SM. 

In this letter, we propose to overcome this problem 
by using B d -> K^°K^° decays to predict S(B S -> 
K*°K*°) within the SM, including the theoretical error 
associated to hadronic uncertanties, in particular to the 
GIM-suppressed penguin contractions. Conside ring the 
experimental precision that is expected at LHCb [1 11 ] and 
at a future super B-factory [12j |. we give a theoretical 
estimate of the deviation from zero of S(B S — * K*°K*°) 
within the SM. This is a crucial ingredient to search for 
NP effects in this decay mode: a deviation from zero 
much larger than the estimated SM error would be a 
strong signal of NP. 

For a given polarization of the final state (26[, we can 
write the decay amplitude of B s — > K*°K*° decays as 



A(B S 



K* U K 



-v; b v ts p s 



V* V pGIM 
v u b y us s 



(2) 



in the same notation of Eq. {T]). Comparing Eq. |2]) to 
Eq. (|TJ) , it is clear that the same NP in b — » s penguins 
enters both this channel and the golden mode Bd —> 4>Ks- 

From an experimental point of view, the K*° 
can be reconstructed as K*° — > K + ir~ and K 
K~tt + . Since the final state is a CP eigenstate, it is pos- 
sible to measure the CP asymmetry parameters S and 
C from the time-dependent study of the tagged decay 
rates. The information on the flavour of the decaying B s 



mesons 
*o 
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is provided by the usual tagging techniques. Since in this 
case the B s meson decays only to charged tracks directly 
originating from the vertex, the reconstruction and ver- 
texing of the B s mesons should be possible at LHCb, al- 
lowing to measure the parameters of the time-dependent 
CP asymmetry. 

With the same approximation of B d — * <t>Ks, i.e. 
neglecting the CKM suppressed contribution of P^} M , 
the SM expectation values for the coefficients of the CP 
asymmetry are simply given by S(B S — ► K*°K*°) = 
and C{B 3 -> K*°K*°) = 0, as 



Xcp(B s - K*"K*°) = e 2 ^ 



s A(B s -> K*°K*°) 
A(B S -► K*°K*°) 



= 1. 



This is a null test of the SM, but an estimate of the error 
induced by neglecting P GIM is needed. 

The advantage of this mode, with respect to the case 
of Bd — » (j)K , is represented by the possibility of cal- 
culating the theoretical error in a model independent 
way, using the measurement of BR(B S — > K*°K*°) and 
C(B S — > K^K* ), together with the information on 
the order of magnitude of P GIM provided by the time- 
dependent study of B d — > K*°K*° decays. 

The idea follows the calculation of the error on sin 2/3 
in Bd — > J/ipK° presented in ref. The expression 

for the decay amplitude of Bd — > K*°K*° in the same 
notation of Eqs. (JTJ) and ([2} is given by 



A(£? d -> iT u iT u ) = -V^VtdPd - V^VudPr M , (3) 

which is equivalent to Eq. @ , except that in this case the 
two combinations of CKM matrix elements have the same 
order of magnitude. As a consequence, the sensitivity to 
P d is maximal in this case. From the measurement 
of the BR and the CP parameters S and C, fixing the 
CKM elements to their SM values obtained by the UT 
fit one can determine \P d \, P GIM |, and the relative 
strong phase 5d- In the SU(3)-symmetric limit, P^ GIM ^ = 



P 



(GIM) 



and Sd 



Imposing these relations, as done 



in ref. [15], would introduce a (difficult to estimate) error 
associated to SU(3) breaking To be conservative, 

we instead allow for a SU(3) breaking up to 100%, much 
larger than any known breaking effect. 

Starting from these considerations, the estimate of 
the theoretical expectation for the deviation of S(B S — > 
K*°K*°) from zero proceeds through three steps, in anal- 
ogy to ref. [13] : i) a fit to determine P s from BR(B S — > 
K*°K*°); ii) a fit of \P d \, |P d GIM |, and 6 d from the ex- 
perimental values of BR, S, and C of B d -> K*°K*°. In 
the fit, only the solution that gives \P d \ compatible with 
P s | is considered; iii) a fit of the B s — > if* ^* decay 
amplitude from the experimental values of BR and C, 
performed forcing the absolute value P GIM | in the range 
obtained allowing 100% SU(3) breaking effects around 
the central value of |P GIM |. To be conservative, no infor- 
mation from Bd — > K*°K*° is used for constraining \P S \ 
and <L. 





channel 


BR 


S 


C 


B a 




(11.8±0.6)1(T 6 


-0.07 ±0.02 


0.01 ±0.02 


B d 




(5.00±0.25)10~ 7 


-0.12 ±0.02 


0.13 ±0.02 



TABLE I: Input values used to estimate the precision on the 
determination of arg(Acp). 



Let us first discuss the experimental prospects. Based 
on these, we then present an example of how our method 
might work once the relevant modes will be measured. 
An estimate of the level of precision reachable at LHCb 
is difficult to give at this stage and goes beyond the pur- 
pose of this work, since details on the reconstruction of 
the LHCb detector are needed. In addition, the lack of 
measurements of Bd — > K*°K*° makes any prediction 
harder. Nevertheless, few educated assumptions might 
help us to understand the order of magnitude of the ex- 
perimental error on S(B S — ► K*°K*°). We assume that 
i) LHCb will provide a measurement of the S and C pa- 
rameters with an error of ~ 0.02 (comparable to what is 
expected for B s — > K + K~); ii) a 5% precision on the de- 
cay rate will be obtained at LHCb or at a super B-factory 
running at the T(5S I ) resonance 17|; iii) a similar pre- 
cision will be available for B d — > K*°K*° rates and CP 
asymmetries. Concerning this last point, it is important 
to stress that B d — > K*°K*° decays are CKM suppressed 
with respect to B s K*°K*°. Nevertheless, with LHCb 
integrating more than two years of data and/or a super 
B-factory integrating > 30 ab _1 , there should be no lim- 
itation given by the available statistics. For the central 
values, we assume that they lie in the ballpark of the 
calculation of ref. [HI], but we have checked that larger 
values of the BR's give similar results. The values we 
assume are summarized in Tab. U [28]. 

We now give an example of the precision we might ex- 
pect on the theoretical prediction of S(B S — > K*°K*°) 
using the numbers given above. Using our method we ob- 
tain the distribution of S(B S -> K*°K*°) shown in Fig.Q] 
with an RMS of 0.015. This corresponds to a theoretical 
error on arg(\ G p(B s -> K*°K* )) of 0.8°. Clearly this 
estimate is only illustrative as it is based on the values 
in Tab. U inspired by factorization models. Once data 
will be available, however, the method will provide an 
estimate independent of any theoretical model. 

The presence of multiple polarizations in the K*K* fi- 
nal state does not change the idea we propose here, but 
it has a practical impact on the analysis strategy. Our 
procedure can be followed for each polarization (longi- 
tudinal or transverse), taking into account the relative 
minus sign in the CP eigenvalue. Experimentally, using 
an angular analysis it is possible to separate the different 
contributions and independently determine rates and CP 
parameters [29| . 

In terms of the experimental fit, the separation of the 
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FIG. 1: Probability density function for S (top) and arg(Acp) 
(bottom) for the decay B s — > 7f*°/("* obtained with the pro- 
cedure detailed in the text, using the input values given in 
Tab. U 



three polarizations requires to add the angular distribu- 
tion of the final state particles to the maximum likelihood 
fit, as it was done for example in the Bd — > p + p~ time- 
dependent analyses of BaBar [2(| ■ 

From a practical point of view, the presence of three 
polarizations helps to increase the experimental preci- 
sion, with respect to the case of a single polarization. 
In fact, the measurement is also sensitive to interfer- 
ence terms among the different polarizations, as for 
Bd — ► J /ipK* time-dependent analyses performed at the 
B-factories [2l], [HI • One can define eleven parameters de- 
scribing the complex Bd decay amplitudes for the three 
polarization states and their CP conjugates, up to an 
arbitrary global phase. In terms of these eleven parame- 
ters one can compute three sets of rates, S and C coef- 
ficients (one for each polarization). Since the number of 
unknowns is smaller than the number of observables, the 
presence of three polarizations in the final state will im- 
prove the precision of the analysis. The same procedure 
can be used for B s decays. 



In principle, the same approach can also be applied to 
B s -> K*°K°, B s -> K*°K° and B s -> K°K° decays, 
with the caveat that the strategy has to change in order 
to face the different experimental challenges. 

For B s — + K°K°, the measurement of the BR should 
be possible at LHCb or at a super B-factory. On the 
other hand, the time-dependent CP parameters <S* and 
C cannot be measured, since the extrapolation of the 
B s vertex from the flight direction of two Ks does not 
seem possible at LHCb, while a B-factory has not enough 
vertex resolution to follow the fast oscillations of B s 
mesons [I?} • Nevertheless, it is still possible to obtain a 
determination of Xcp, measuring the tagged decay rates 
for At > and At < 0. The sign of At can be mea- 
sured at a super B-factory, using the Ks flight direction 
to determine the B vertex [23| • Using a full Monte Carlo 
simulation, it was shown that it is possible to measure 
arg(Ar7 p {B s — > K°K )) with an experimental error less 
than 20° [17|. The actual error could be even smaller, 
if the improvement of the vertexing detector (due to the 
use of a layer zero of the silicon detector close to the 
beam pipe) will allow to separate primary and secondary 
vertices on B — » DX decay s 1241 . strongly reducing the 
background contamination [301 ] - At the same time, the 
RMS for arg(Acp) expected in the SM, taking into ac- 
count P S GIM with our method, should be at the level of 
4° 0. 

The case of B s -v K_*°K° and B s -> K*°K° is more 
similar to B s — > K*°K*°. The main difference in this 
case is that there are two different particles in the final 
state. As a consequence, the number of hadronic pa- 
rameters to determine is twice the number of hadronic 
parameters for a single polarization in B s — > K*°K*° 
modes. On the other hand, the number of experimental 
observables is larger. Reconstructing B s — > K*°K Q and 
B s -> K*°K° (B„ -> K°K*° and B s -> K°K*°) from 
K + ir~Ks (K~ir + K$) final states it is possible to mea- 
sure CP violating effects [§5j , which provide four observ- 
ables (S, C, S, and C) in addition to two decay rates. It 
will be possible to use B s -> K*°K° and B s -> K°K*° to 
obtain two null tests of the SM, using the upper values on 
the two P S GIM contributions obtained from Bd —> K*°K° 
and B d — > K°K*° respectively. 

To summarize, we have proposed a new strategy to 
look for NP in b — ► s penguins without relying on model- 
dependent estimates of the hadronic uncertainties. The 
new golden channel we suggest is B s — > if(*)°J^(*)°. 
We claim that the SM pollution in the null tests of 
the SM from time-dependent CP asymmetries in this 
golden channel can be controlled with a high accuracy 
in a model-independent way. The key observation is 
that, even allowing for SU(3) breaking effects of 0(1), 
using the experimental information on the SU(3)-related 
channel B d K^K^ it is possible to put a strong 
constraint on the polluting CKM-suppressed penguin 
amplitude. The most promising channel seems to be 
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B s — > if* ^* , which can be reconstructed from four 
charged tracks in the final state and should be easily 
accessible at LHCb, together with the B d decay to 
the same final state (which can already be studied 
with the full dataset collected by BaBar and Belle). 
Pseudoscalar-vector and pseudoscalar-pseudoscalar final 
states imply the presence of K$ mesons, making the 
analysis harder in the environment of a hadron collider. 
In this respect, a super B-factory would play a very 
important role. 

We warmly thank Y. Xie for pointing out to us an 
embarassing error in the previous version of this paper. 
Fortunately correcting this error further strengthens our 
proposal. 
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of the final state. For simplicity, we omit the correspond- 
ing label. We comment on the experimental aspects of 
measuring vector final states and of extracting polarized 
amplitudes after presenting our method for a single po- 
larization. 

In ref. ph| ]. SU(3) was used in conjunction with QCD 
factorization to relate CP violation in B s and Bd de- 
cays to two kaons. This is another interesting example of 
model-dependent predictions for CP violation in b — ► s 
penguins. 

During the review process of our manuscript, the BaBar 
collaboration published a first measurement of BR(Bd — * 
K*°K*°) = (4.9±i;| ± 0.5) ■ 10" 7 , in agreement with the 
value we assumed, altough with a large uncertainty [TjJ. 
In case large correlations are observed among the three 
rates and the three sets of CP parameters, the theoretical 
analysis will have to take them into account. 
Concerning the tagging at the T(5S), the training of a 
tagging algorithm for the B s will face new challenges, as 
both K and K can be produced in a B 3 decay. One could 
cope with this problem using, for instance, the difference 
in the momentum of the kaons coming from the b quark 
decay and of the kaons coming from the cascade of D 3 
mesons produced by the B s . A detailed estimation of the 
change in performances goes beyond the purpose of this 
paper. We use the current B-factory performance as a 
crude estimate of the tagging efficiency. 



